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Abstract 

The problem of controlling and stabilising solutions to the Kuramoto-Sivashinsky 
equation is studied in this paper. We consider a generalised form of the equation in 
which the effects of an electric field and dispersion are included. Both the feedback 
and optimal control problems are studied. We prove that we can control arbitrary 
nontrivial steady states of the Kuramoto-Sivashinsky equation, including travelling 
wave solutions, using a finite number of point actuators. The number of point actua¬ 
tors needed is related to the number of unstable modes of the equation. Furthermore, 
the proposed control methodology is shown to be robust with respect to changing 
the parameters in the equation, e.g. the viscosity coefficient or the intensity of the 
electric field. We also study the problem of controlling solutions of coupled sys¬ 
tems of Kuramoto-Sivashinsky equations. Possible applications to controlling thin 
film flows are discussed. Our rigorous results are supported by extensive numerical 
simulations. 


1 Introduction 

The Kuramoto-Sivashinsky (KS) equation on L —periodic domains 

Ut T u X xxx T Uxx T uu x — 0, (Ll) 

u(x, t ) = u{x + L, t), 

is a paradigm evolution equation that has received considerable attention in recent years 
due to its wide applicability as well as the rich and complex dynamics that it supports. 
The KS equation arises in many physical problems including falling film flows [7, 62, 59, 
29], two-fluid core-annular flows [50, 19], flame front instabilities and reaction-diffusion- 
combustion dynamics [60, 61], propagation of concentration waves in chemical physics 
applications [42, 43, 41], and trapped ion mode dynamics in plasma physics, [15]. The 
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KS equation (1.1) is one of the simplest partial differential equations (PDEs) that can 
produce complex dynamics including chaos - see for example the numerical experiments 
in [30, 31, 33, 39, 52, 63, 74, 75]. Routes to chaos have been shown numerically to follow 
a Feigenbaum period-doubling cascade - see [63] where the two universal Feigenbaum 
constants are also computed for the KS with three-digit accuracy. A detailed knowledge 
of the stationary, travelling and time-oscillatory solutions (typically chaotic) of (1.1) is 
significant in technological applications that seek to enhance heat or mass transfer, for 
example. In this sense certain solutions are better than others and a description of the 
solution phase space is a crucial step in constructing relevant control strategies that can 
access unstable states, for instance, that may be desirable in applications. 

In many studies equation (1.1) is scaled to 27T—periodic domains according to the 
rescaling 


* 2?r ,* ( 2 ^\ L 2tt * 2tt 

x = — x ’ f u = 5 V = ( L2 ) 


to take the form (we drop the stars and use the same symbols for dependent and inde¬ 
pendent variables) 


Ut A vuxxxx 4“ v*xx 4~ mix — (1-3) 

u(x, t ) = u(x 4- 27t, t), 

where v — (2 tt/L) 2 is a positive parameter that decreases as the system size L increases. 
The mathematical interest in the KS equation (and related models - see below) resides in 
the fact that it is a simple, one-dimensional equation exhibiting complex dynamics mak¬ 
ing it amenable to analysis and also a good case study in the area of infinite-dimensional 
dynamical systems and their control. The equation is of the active-dissipative type 
and instabilities are present depending on the value of v. If v > 1, it is well known 
[56, 58, 64, 65] that the zero solution representing a flat him, is unique. However, when 
v < 1 the zero solution is linearly unstable and bifurcates into nonlinear states includ¬ 
ing steady states, travelling waves and solutions exhibiting spatiotemporal chaos - the 
dynamical complexity increasing as v decreases. Some of these solutions are stable, and 
others are unstable [39]. In [24, 39, 51], one can find studies of the stability of steady 
states of the KS equation. 

There is an extensive literature on the behaviour of the solutions to the KS equation. 
Well posedness of solutions is studied, for instance, in [56, 64, 65]. It was proved in [18] 
that the long time dynamics of the KS equation are finite dimensional in the sense that 
they are governed by a dynamical system of finite dimension which is at least as large 
as the number of linearly unstable modes (this number scales with L or z/ -1 / 2 for (1.1) 
or (1.3), respectively); these authors also proved that the solutions are attracted by a 
global attractor, a set of finite dimension. Boundedness of solutions for general initial 
conditions was proved independently and by using distinct methods by [17, 27, 32], 
These studies also focussed on finding bounds for the dimension of the global attractor 
by estimating L 2 —norms of the solutions, starting with the odd-parity results of [46] and 
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those for general initial data by [17] along with more recent improvements in [9] and [47]. 
Analyticity of solutions in a strip in the complex plane around the real axis was also 
proved in [16] and [3] using different methods. 

In the context of falling him hows there have been several studies to extend the KS 
equation by including additional physical effects. Of most interest to the present study 
are the derivations in [70, 72] for him how over hat walls in the presence of electric helds 
applied perpendicular to the undisturbed interface. The resulting equation, that also 
incorporates the effects of dispersion, is a generalisation of (1.3) and takes the form 

Uf + VUxxxx + liH[ ^XXX ] T Sttxxx T uxx T uu x — 0, (1-4) 

u(x, t) = u{x + 2-7T, t), u(x, 0) = Uo(x), 

where p > 0 measures the strength of the applied electric held and the parameter 5 
measures dispersive effects. The linear operator H is the Hilbert transform operator and 
represents how destabilisation due to the electric held. On 27T—periodic domains the 
definition of H is 

n[u\(x) = PV u(0 cot ( L5 ) 

where PV stands for the Cauchy principal value integral. In the model analysed, the elec¬ 
tric held needs to be found by solving a harmonic problem above the him and calculating 
the Dirichlet to Neumann map of the solution to construct the Maxwell stresses that 
interact with the hydrodynamics - see [70, 72] for the details. In fact, for the linearised 
problem the eigenvalues A corresponding to the eigenfunctions exp (ikx) are 

A = k 2 + /ik 2 \k\ — uk A + i8k 3 , (1.6) 

showing that the presence of the electric held destabilises the how and increases the num¬ 
ber of linearly unstable modes. Note that instability is possible if \k\ < k c = +4 -, 

and so there are 21 + 1 unstable modes where l is the integer part of k c (this fact is used 
in Proposition 1 later). This additional destabilisation is important in what follows and 
makes the control problem more challenging. The modified equation (1.4) in the absence 
of dispersion (<5 = 0), has a similar dynamical behaviour to the KS equation (1.3) but 
with chaotic dynamics appearing at higher values of v as fi increases. In fact, bound¬ 
edness of solutions and an estimate of the dimension of the global attractor have been 
proved in [71] for a class of more general operators whose symbols in Fourier space are 
such that the electric field term in (1.6) is \k\ a with 3 < a < 4. On the other hand, in the 
absence of an electric held but with dispersion present, it is established that dispersion 
acts to regularise the dynamics (even chaotic ones) into nonlinear travelling wave pulses 
- see [36, 37, 2] as well as [53, 69] for a weak interaction theory between pulses that are 
sufficiently separated. 

More recently, Christohdes et al. [4, 6, 10, 11, 12, 21, 45] showed how to stabilise 
the zero solution of the Kuramoto-Sivashinsky equation by using state feedback controls. 
They also proved that using linear feedback controls, it is possible to stabilise the zero 
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solution using only 5 point actuated controls. In addition they prove that the stabil¬ 
isation is possible using only a certain number of observations of the solution instead 
of full knowledge of the solution at all times, as long as the number of observations is 
equal to or exceeds the number of unstable modes. In further work utilising nonlinear 
feedback controls [4, 11], Christohdes and co-workers formulated optimisation techniques 
and computed possible optimal states by analysing a large number of runs; a proof of 
the existence of these optimal positions was not given, however. 

In this work we use linear feedback controls and techniques similar to those in [6, 11], 
to stabilise non-uniform unstable steady states of generalised versions of the KS equation 
at small values of v that have not been attempted yet. The mathematical complication 
is due to the increase of the number of linearly unstable modes as v decreases and /i 
increases - see (1.6). We achieve stabilisation of non-uniform states by stabilising the 
zero solution of a modified PDE that is satisfied by a perturbation to the desired steady 
state. The resulting equation to be controlled becomes 

rri 

Ut + VU X xxx T T 3'U‘xxx T u xx + uu x — ^ ^ bi(x)fi(t'). 

i= 1 

All our results are still valid in this as well as other cases, as long as the linear operator 
of the PDE has a self adjoint part, a well defined separation between stable and unstable 
modes, and a bounded nonlinearity Af(u) in an appropriate functional space. Such 
conditions are fairly generic in physically derived systems and pose little restriction to 
our methodology. 

In applications there may be some uncertainties in the estimation of the parameters of 
the equation, for example, if the intensity of the electric held or the dispersion parameters 
are not known exactly. It is important, therefore, that the controls applied are robust, 
that is they still work even when these uncertainties are present. We use results from 
control theory [35] to prove analytically that the controls are robust to uncertainties in 
v, 6 and /x, as long as the error in the prediction of the parameters is small enough, and 
present numerical simulations that demonstrate this point. 

A natural question to address after robust stabilisation of the zero solution to the 
modified PDE is achieved, is whether this can be done in an optimal manner. By this we 
mean stabilisation while minimising a cost functional that measures how close we are to 
the desired solution and how much energy we are spending with the controls. This cost 
functional is of the form 

rp rp fJX 

c ( u i F ) = \[ ||w(-,t) - u \\ 2 dt + h\u(-,T) -u|| 2 +! f ^ ^fi{t) 2 dt , (1.7) 

Jo Jo i=1 

where T is the final time of integration and u is the desired steady-state we are con¬ 
trolling. As our control variables we will consider the positions of the control actuators, 
following Lou and Christohdes [45]. Note that the presence of a nonlinearity in the PDE 
makes the problem non-convex and therefore we do not expect an optimal control to be 
unique. However, we only wish to prove the existence of an optimal control and to find 
computationally such optimal controls. 
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The methodology developed and implemented here can also be applied to systems 
of nonlinear coupled PDEs. Of particular interest are systems of coupled Kuramoto- 
Sivashinsky equations that arise in applications to interfacial fluid dynamics problems. 
Such equations were derived systematically using asymptotic methods in [49] to model 
the nonlinear stability of three immiscible viscous fluids of different properties flowing in 
a stratified arrangement in a plane channel under the action of gravity and/or a driving 
pressure gradient. The ensuing dynamics is very rich and in fact instabilities can emerge 
even in the absence of inertia, unlike analogous two-fluid flows. Coupled nonlinear sys¬ 
tems are mathematically significantly more challenging than scalar PDEs since analytical 
results on global existence and estimates of solution norms, for example, are poorly un¬ 
derstood. Detailed computational results into the complexity of the solutions (especially 
their zero diffusion limits) of such coupled systems of KS equations can be found in [48]. 
In the present study we consider the problem of feedback and of optimal control when 
the equations are coupled through the second derivatives alone. This is a special case but 
arises in the derivation of the equations, see [49]; more generally the nonlinear terms are 
also coupled and can cause hyperbolic-elliptic transitions by supporting complex eigen¬ 
values of the nonlinear flux functions, see [48]. We find that solutions to such a system of 
ocupled PDEs can also be controlled through a linear feedback loop and also optimally. 






Figure 1: Bifurcation diagram of the L 2 -norm of the steady state solutions (solid curves 
- blue online) and travelling wave solutions (dashed curves - red dashed online) to the 
gKS equation (1.4) (with 5 = 0) in the presence of an electric field for 0.1 < v < 1, and 
fi — 0 (la), fi — 0.2 (lb), fi = 0.5 (lc), and p, = l(ld). 

The rest of the paper is organised as follows. In Section 2 we prove rigorously that 
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nontrivial solutions to the generalised Kuramoto-Sivashinsky equation can be stabilised 
and analyse the robustness of the proposed controls to variations in the parameters of 
the equation and in Section 3 we present numerical experiments that confirm our results. 
In Section 4 we prove that there exist optimal distributed controls for the KS equation, 
construct an algorithm to optimise the placement of the control actuators and show our 
numerical results. Finally, in Section 5 we extend our results to a system of coupled KS 
equations. We discuss our results in Section 6. 

2 Stabilisation of nontrivial unstable steady states using linear feedback 
controls 

Our goal is to stabilise non-uniform unstable steady states or steady-state travelling 
wave solutions of equation (1.4). For the theoretical analysis of the feedback control 
problem for the KS equation we need L°° bounds on the solution and its derivatives. To 
establish such estimates we will use well known L 2 bounds, together with the Sobolev 
embedding theorem. Optimal estimates for the solution of the KS equation (1.1) in (0, L) 
were obtained by Otto [47], and for the re-scaled 2n —periodic KS equation (1.3) these 
estimates can be expressed in terms of u — (qy) 2 to find 


limsup ||u(-,t)|| < 0(u 1 / 6 ), 

(2.1a) 

t—>QG 


limsup \\u x (-,t) || < C> (V/ 2 ln 5 / 3 (V 1 / 2 )) , 

(2.1b) 

t^oo \ \ J J 


limsup \\u xx (-,t) || < O (is In 5 / 3 , 

(2.1c) 


t—>oo 


where || • || = dx ) 1 / 2 denotes the L 2 -norm of the solution. For the generalised 

equation (1.4) but in the absence of dispersion (<5 = 0), Tseluiko & Papageorgiou [71] 
use the backround flow method to obtain similar estimates in the presence of an electric 
held. Their estimates are of the form 

IM| < (||wo|| + |M|) e Dt + C(u, p) + ||(p||, (2.2) 

where C and D are constants depending on /i and v. and <p is a constructed function with 
finite L 2 -norm (we do not need to give it here). They also proved that for uq € H^{ 0, 27t), 
the first and second derivatives of the solution are bounded, and therefore, u G i7 2 (0, 2n), 
where H* is the Sobolev space of s— times differentiable functions that are periodic and 
have zero mean. 

It is also important to remark that in the case of the generalised KS equation (1.4), 
with both S and /i non-zero, it was proved in [22, 23] that the L 2 norm of the solution is 
also bounded. In fact 


limsup ||u(-, t)\\ < 

t —>oo 


0(v-' 1 !*), if v< I/O, 
C, if v > uq, 


(2.3) 
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where pq depends on the symbol of the linear operator. Note that these are not optimal 
bounds. The authors also prove boundedness in L 2 of spatial derivatives of u up to order 
4. The estimates (2.1)-(2.3) imply (by use of the Sobolev embedding theorem) that there 
exist constants Ci, C 2 that depend only on v and p such that 

Halloo T C , i||'u||j^2, llw^lloo < C 2 1| || H 1 • (2.4) 


The controlled generalised KS equation can now be introduced and will form the basis 
of our analysis and computations. This consists of a forced version of (1.4) and reads 


Uf T VUxxxx T I4'ht\u x xx\ T 3u X xx T u xx T uu x — ^ ^ bi(x)fi(t), x G (0, 27f), t > 0, 

i=i 

x G (0, 27t), 

7/ 

.(or. 4 I- 27 r. i\ = - ( or. tV or. f 1 (C 

dxJ 


u(x, 0) = Uq{x), 

7 / 7 / 

-(x + 2t T,t) = Q-j(x,t), x G (0,2tt), t > 0, 

fi(t) 6 T(0,T). 

(2.5) 

We assume that the initial condition satisfies uq G Hp(0,2n), m denotes the number of 
controls, bi(x), i — 1 ,,m are the control actuator functions and fi(t), i — 1 ,,m 
are the controls. We will use point actuator functions, which means that the functions 
bi{x) are delta functions centered at positions ap, i.e. bi{x) — 5(x — X{), or a smooth 
approximation of such delta functions. 

We use an argument similar to [5, 10, 12] to prove that it is possible to stabilise 
nontrivial steady states of the generalised KS equation (1.4). Using the Galerkin repre¬ 
sentation of u, 

OO OO 

u(x, t ) = u s n (t) sin (nx) + w^(t) cos( 


n=l 


n =0 


nx), (2.6) 

substituting into (2.5), and taking the inner product with the functions ^==, sll ^U r ) and 
e°s(nx) ^ n _ ^ ^ oo, we obtain the following infinite system of ODEs (dots denote time 

derivatives): 


= (— vn A + pn 3 + n 2 ) u s n + 8n 3 u c n + g s n + ^ b s in fi(t) n — l,...,oo 


2=1 

m 


(2.7) 


= (— vn A + pn 3 + n 2 ) u r n — dn 3 u s n + g r n + ^ K n fi(t) n — 0,... ,00 


i=i 


where bf n — J 0 2?r bi(x) sm(nx)dx and b c in — n b{(x) cos(nx)dx. The nonlinearities g( 
and g £ are given by [1] 

9n ^ y ('U^'U^, ^j^k) + 2 ^ y ('Uj'U^, + ^ ‘ ‘ ‘ ’ 


-27T 


j+fc=n 

n 


j—k=n 


uC j u k + ~ w j?4), 71 = °, • • •, oo- 

j+k=n j—k=n 


2y/7V . 
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In deriving the system (2.7) we used 77[sin(x)](x) = — cos(x) and H[ cos(x)](x) = sin(x); 
these formulas can be derived using contour integrations in the complex plane, for ex¬ 
ample. 

We now define 


£ 


U 


5 


where 


Zu = K uf u\ ■ ■ ■ uf uf] T , z u s = [uf +1 


H+i 


where z™ contains the coefhcients of the (slow) unstable modes and zf those of the (fast) 
stable modes. In addition, 


G = [0 g{ g\ g s 2 g\ ■ ■ ■ } T , D = diag (0, 5, -5, ■ ■ ■ , Sn 3 , -Sn 3 , • • •) , F = [hit) f 2 (t ) 
Furthermore, we introduce the notation 


A = 


A u 

0 


0 

A s 


and B — 


B s 


( 2 . 8 ) 


where 


A u = diag(0, -v + n + l, -v + n + l, ■ ■ ■ , -l A u + /x/ 3 + l 2 , -l A v + /x/ 3 + l 2 ), 

A s — diag(—(/ + 1 ) 4 zx + /x(Z + l) 3 + (/ + l) 2 , — {l + 1 ) 4 zx + /x(/ + l) 3 + (/ + l) 2 , • • •), 



r h c 
°10 

h c 

°20 

... h c 

y m0 





h s 

°11 

h s 

°21 

... h s 

y ml 


r b s 

b s 

• b s 1 





. . . b c 


°i(/+i) 

°2(/+l) 

U m(l+ 1) 


B u — 

u n 

U 21 

u ml 

, B s = 

h c 

h c 

°2(Z+1) 

°m(l+ 1) 

• (2-9) 


h s 

°ll 

h c 

°2l 

... h c 

u ml 


- 


' 



h s 

L °u 

h s 
°21 

... h s 

U ml J 





We can rewrite the infinite dimensional system of ODEs (2.7) as 

z u = Az u + Dz u + G + BF. (2.10) 


We have the following result. 

Proposition 1. Let u be a linearly unstable steady state or travelling wave solution of 
(1.4) and let 21 + 1 be the number of unstable eigenvalues of the system 

Uf — VUxxxx lAH-\u X xx\ 'Uxxi ( 2 . 11 ) 

i.e., I + 1 > M+V 2 ^, +4 ~ > L If m = 21 + 1 and there exists a matrix K £ R mXm suc h 
that all of the eigenvalues of the matrix A u + B U K have negative real part, then the state 
feedback controls 

l/l ••• fmf=F = K(zl-zl) (2.12) 

stabilise u. 
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Proof. Let u — u + v be a solution to (1.4). Substituting into (1.4) and using the fact 
that u is a steady state or travelling wave solution of (2.5), we obtain the following PDE 
for the perturbation v, 

Vt + w xxxx 4“ pB[v XX x\ 4“ $v xxx 4~ v xx 4~ vv x 4~ (uv) x — 0; (2.13) 

and in controlled form we have 

m 

vt 4~ vv XX xx 4~ f^B[v xxx \ 4~ 8v xxx 4~ v xx 4~ vv x 4~ (uv) x — ^ ^ bi(x)fi(t). (2.14) 

i=l 

First we will prove that the given controls stabilise the zero solution of 

Vt — vv xxxx 4~ lAH.\v xxx \ 4~ v xx . (2.15) 


Note that the dispersion term does not affect instability and it is not necessary to include 
it in this part of the analysis. After applying a Galerkin truncation and the controls given 
by (2.12), we obtain 

■ v _ A u 4- B U K 0 
2 “ L B,K a s 

Since the eigenvalues of A u + B U K have negative real part and the matrix multiplying 
z v is triangular, it follows that the zero solution to (2.16) is exponentially stable. 

Next, following [5, 10, 12], we use a Lyapunov argument to show that these controls 
stabilise the zero solution to equation (2.13). We first use the fact that exponential 
stability of the system (2.16) implies that there exists a positive constant a such that 
the operator A{v) = -vv xxxx ~ lM[v xxx ] - v xx + Ya =i h(x)K l .zf, where K v denotes the 
i— th row of the matrix K, satisfies 


= Cz l 


(2.16) 


(Av, v) < —a||u|| 2 . 


(2.17) 


Defining V(v) = J Q 27r ^ dx, it is easy to verify that V(0) = 0 and V ( v ) > 0, Vu > 0. 
Multiplying (2.13) by v and integrating gives 


d f 2n v 2 , 

— / — dx — 

dt 


-27T 


-27T 


-27T 


-27r 


'0 


vv t dx = (Av, v) — 8 v xxx v dx — v 2 v x dx — v(uv) x dx. 

Jo Jo Jo 

(2.18) 

Integration by parts and and use of periodicity shows that the first two integrals on the 
right-hand side of (2.18) are zero. It remains to obtain an estimate for the third integral. 
Again using integration by parts and periodicity gives 


-2tt 


-2tt 


-2tt 


(uv) x v dx — — uv x v dx — 

Jo 

inf u 3 
< -- 


2 ; 1 

u r v dx = — 


-27r 


u T v 2 dx 


-2tt 


inf u 3 


Ikll 2 - 


2 


v 2 dx = 


2 



Adding everything np, we obtain 


1 d 

2 dt 



< - 




(2.19) 


If the eigenvalues of the matrix A u + B U K are chosen such that 2 a + inf u x > 0, we 
obtain that ^V{v{t)) < 0, which proves that V is a Lyapunov function for the system 
and therefore the zero solution is stable. 

Using controls (2.12), we can therefore stabilise the nontrivial steady state u of the 
original equation. □ 

Using Proposition 1, we can conclude that in order to stabilise the steady state u of 
equation (1.4) we should solve the PDE 


m 

Ut 4“ nu xxxx 4“ lBht[u X xx\ 4~ $'U j xxx 4~ w xx 4~ uux — ^ ^ bi(x')K{.(z u Zy). (2.20) 

i =1 

Remark 1. Since the solutions to the generalised Kuramoto-Sivashinsky equation are 
taken to be periodic with mean zero, it follows that inf u x < 0. Therefore, the constant a 
in (2.19) must be chosen large enough so that a 4- mi ™ x is positive. In the case when 6 > 0, 
we also need to account for the fact that the amplitude of the solutions (and therefore the 
absolute value of their derivatives) grows with 6 [38, 2j. Further details can be found in 
Section 3. 


Remark 2. The above proposition is clearly valid for the case when u — 0, in which 
case the controls are fi(t) — K^zf, as presented in [5, 10, 12J. In the case when u is 
a travelling wave, the result follows using a time dependent z u . See Section 3 and in 
particular Equation (3.6). 

Remark 3. From estimates (2.1) ; it follows that the value inf \u x \ is finite and therefore 
we can conclude (2.19). 

Remark 4. Christofides et al. [5, 10, 12j argued that due to the multiplicity of the 
eigenvalues of the Kuramoto-Sivashinsky equation being less or equal than 4, one would 
only need 5 controls to stabilise the zero solution of that equation. The same holds in our 
case: the multiplicity of the eigenvalues of the linear operator in (2.15) is also less than 
or equal than f, but numerical results suggest that we need to use m = 21 4- 1 controls, 
or at least a number of controls that is close to the number of unstable modes, see Fig. 2 
and the discussion below. 


Remark 5. The fact that we are separating the system between stable and unstable modes 
implies that the matrix B u is square (B u e W mXm ), and using bi{x) — S(x — Xi) means 
that B u has full rank. It follows that the Kalman rank condition [76] is automatically 
verified and the matrix K needed for the stabilisation will always exist. 
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2.1 Robustness of controls 


A natural and important question is whether the proposed control methodology is robust 
with respect to changes (or uncertainty) in the parameters v. p and 5 that appear in the 
equation. The robustness of our method can be proved rigorously using techniques from 
control theory, e.g. [35, Thm. 6], and we take this up next. 


Proposition 2. Let Aq i — be the eigenvalues of the matrix C appearing in 

(2.16), X be the matrix of eigenvectors of C , and let k(-) denote its condition number. 
Then we have 


K lb < 


( ll^lh + max (|Aj|) k(X) \ 


&m(B) 


( 2 . 21 ) 


where a m (B) is the m-th smallest singular value of B, which is defined in Equations (2.8) 
and (2.9), and the solution z v to equation (2.16) satisfies 


|^ v (t)|| < ft(X)max 
3 


A it 


'( 0 )||. 


( 2 . 22 ) 


Proposition 3. If the feedback matrix K is such that Equation (2.16) is exponentially 
stable, then the perturbed closed loop system matrix A + BK + A remains stable for all 
disturbances A which satisfy 


A ||2 < mimjyv (si — (A + BK)) —: ((K), (2.23) 

s=iuj 


where 

C(A’) > min Re (^) . 
and Re(-) denotes the real part. 

In particular, if there is an error in the estimation of the parameters v and p, then 
the feedback matrix K will still stabilise the zero solution as long as the error in the 
parameter estimation is bounded by ((K). We have studied the robustness of the controls 
for stabilising steady states and travelling waves by combining Propositions 2 and 3. We 
now present a summary of our results. 

Variations in 5 

As seen from the dispersion relation (1.6), variations in <5 do not affect the stability of 
the solutions and consequently so they do not affect the matrix K. This implies that the 
matrix A is zero, and the zero solution to system (2.16) is still stable. In the case when 
we are interested in stabilising travelling waves, we need to take into account the fact 
that the amplitude of the travelling waves increases with 5 - see for example [38, Fig. 1], 
Hence, if we overestimate the value of inf u x and take this into account when choosing 
the new eigenvalues, the stabilised solution should remain close to the desired travelling 
wave as demonstrated by our numerical experiments in Figure 4. 
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Variations in v and p. 

As seen from (1.6) variations in v and fi can affect the stability of the solutions and the 
number of unstable modes. An increase in unstable modes in turn affects the number 
of controls needed since our theoretical results support that we need the same number 
of controls as unstable modes as stated in Remark 4. However, we have performed 
numerical experiments (see Fig. 2) that show that using two less controls than predicted 
theoretically does not affect the stability of the solutions. 



(a) t — 1 



(c) t = 5 


0.5 

0 



0.5 

0 


j— 0.5 

* 1 



(b) t — 3 



(d) t — 7 




Figure 2: Snapshots of the time evolution of the stabilised travelling wave solution in 
figure 10b using m = 19 controls instead of m = 21 at different times. Red full line is 
the controlled solution, black dashed lines are the desired travelling wave and the dots 
represent the controls and their intensity. 


Now we consider the case where we have some uncertainty of amplitude ei and 62 in 
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the values of v and fi, respectively, 

m 

Ut = -(v + ei) ^XXXX (h + e 2 )'H[u xxx\ UXX - uu x + E bi(x)K t . (z“ - z u ). (2.24) 

i=1 


The controls have been chosen so that the solution to the equation is stabilised when 
ei — 62 — 0. Multiplying (2.24) by u, integrating by parts and using Young’s inequality, 
we find 

“IN -,£)|| 2 < -k|M | 2 - ei\\u xx \\ 2 + y (||u x || 2 + \\u xx \\ 2 ) , 

where k — a + is a constant. On the other hand, the perturbation —e\u xxxx — 

£2 'H[uxxx\ can be discretised and written as 

A = diag(0, — ei/c 4 + 62 k 3 , —e\k A + 62 k 3 ), (2.25) 

k = 1,..., N/2, and it follows that its Frobenius norm is given by 

N /2 N /2 

|| A||2 = 2 k 6 (—ei k + 62) 2 = 2 k 6 (elk 2 — 2ei62 k + e|) . (2.26) 

k=1 k =1 


For stability we need (2.26) to satisfy estimate (2.23) - see Proposition 3. Therefore, we 
have the following proposition. 

Proposition 4. Let K be a matrix such that A U + B U K has the prescribed (negative real 
part) eigenvalues, Ai,..., X m , with m = 21 + 1, and let 


BI< 


B U K 0 
B S K 0 


Then the perturbed system A + BK + A ; where A is given by (2.25), is stable provided 
that 


N/2 


1/2 


2 k 6 (elk 2 — 2ei62 k + 62) < min ctn (si — (A + BK)). 


k =1 


We have performed numerical experiments to test the robustness of the controls, 
and in particular we focussed on robustness with respect to the parameters <5 and v. 
Numerical results are presented in Figures 2-4 (results in these figures as well as in 
Figure 10, are shown in the original unsealed domain of length L - see (1.2) for the 
transformations). In figure 2 we use the same parameter values as in Figure 10b but 
we use 19 controls instead of 21, i.e. two controls less than the number of unstable 
eigenvalues. The dashed curve is the desired travelling wave solution and the solid curve 
(red online) is the controlled solution with 19 controls. We conclude, therefore, that 
our control methodology is robust with respect to a slight decrease in the number of 
controls. Note however, that the number of controls cannot be significantly smaller than 
the number of unstable eigenvalues - for example, running the same numerical experiment 
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with 17 controls did not yield satisfactory results in the sense that wavy perturbations 
observed in panels (b) and (c) were not suppressed. 

A robustness test with respect to changes in v (with <5 = fi — 0) is depicted in 
Figure 3. We begin with an unstable travelling wave at v — 0.013 and wish to control 
it but by solving the KS equation with a reduced value of v = 0.01, i.e. we impose an 
uncertainty in the value of the parameter v or equivalently in the shape of the desired 
solution. The results again show robust behaviour with the two solutions being almost 
indistinguishable. Finally in Figure 4 we present robustness experiments for /r = 0, 
v — 0.01 and changes in the dispersion parameter 8 from 0.03 to 0.04, with equally 
accurate performance as before. 








Figure 3: Snapshots of the time evolution of a stabilised travelling wave solution for 
8 = /i = 0 and assuming uncertainty in the parameter v. Black dashed line is the desired 
travelling wave (which is the correct solution for v = 0.013 -O- L 55) and red full line 
is the controlled solution assuming v = 0.01 L « 62. 
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Figure 4: Snapshots of the time evolution of a stabilised travelling wave solution for 
fi = 0, v — 0.01 and assuming uncertainty in the parameter 5. Black dashed line is the 
desired travelling wave (which is the correct solution for <5 = 0.03) and red full line is the 
controlled solution assuming 5 = 0.04. 


3 Numerical Results 

Section 2 was devoted to proving rigorously that steady states and steady-state travelling 
wave solutions of the generalised KS equation can be stabilised using linear feedback 
controls. The number of controls is predicted to be at least as large as the number of 
linearly unstable modes, and robustness with respect to uncertainty in the parameters 
v and 8 was also proved. In this section we implement the linear feedback controls 
numerically and undertake an extensive computational study of the stabilisation and 
control in practical situations. 
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3.1 Construction of non-uniform steady states and travelling waves 

One of the main objectives of the present work is the stabilisation of unstable solutions 
of the gKS equation. To obtain steady state solutions u{x) (in the absence of dispersion), 
we need to solve the equation 

VUxxxx + /iH[u xxx] T U xx + wu x = 0, (3.1) 

in the interval [0, 27r], subject to periodic boundary conditions. Travelling waves of speed 
c are found by looking for solutions of the form u(x,t) — U{x — ct) — U(^) and solving 

- cU' + vU"" + fiH\U"'] + 8U'" + U" + UU' = 0, (3.2) 

subject to periodic boundary conditions, where primes denote differentiation with respect 
to £. We note that equation (3.1) is a particular case of (3.2). Expressing the solutions 
in Fourier series 

OO 

U(£) = ^E/£sin(n£) + E/£cos(n£), (3.3) 

n=1 

and substituting into (3.1) and (3.2) we obtain an infinite system of nonlinear algebraic 
equations for the coefficients U*, U, n = 1,..., oo, or for the coefficients and the velocity 
c, in the case of travelling waves. The resulting system of equations for steady states is 


(un 4 — fin 3 — n 2 )U.— 0, n — 1,..., oo, (3.4a) 

(un 4 — fin 3 — n 2 )U,£ + fj s n = 0, n — 1,..., oo. (3.4b) 

For travelling waves we can assume, without loss of generality due to translation invari¬ 
ance, that U( = 0, to obtain 

— (cn + 8n 3 )U.n + (un 4 — fin 3 — n 2 )U^ + g c n = 0, n = 1,..., oo, (3.5a) 

(cn + 8n 3 )Un + {nn 4 — fin 3 — n 2 )U.£ + g s n — 0, n — 2,..., oo, (3.5b) 

{c + 8)U{ + g{ = 0, (3.5c) 


The systems were truncated and solved using a nonlinear solver (e.g. Matlab’s fsolve ) 
to find solutions to system (3.4) by first setting fi — 0 and carrying out a numerical 
continuation on u, and secondly by fixing the desired value of v and varying fi. For 
travelling waves we used continuation on v, fi and 8. Without loss os generality we also 
impose c > 0: if U{x — ct) is a solution of (3.2) with c < 0, then —U(—x — (—c)t) is also 
a solution with c > 0. 

Given the Fourier coefficients and the velocity of a travelling wave, we can write the 
soluton of the KS equation as 

u(x,t) = U(x — ct) — (Pn cos(nct) + sin(nct)) sin(nx)+ , . 

E^°=i {Un cos(nct) - U s n sin(ncf)) cos(nx). 

Our computational results are presented in the bifurcation diagram in Figure 1 that 
depicts the variation of the T 2 -nornr with v of the steady states and travelling wave 
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solutions of the gKS equation (1.4) in the absence of dispersion (<5 = 0). Panels (a)-(d) 
correspond to fi — 0, 0.2, 0.5, 1.0; steady-states are plotted with solid curves (blue on¬ 
line) and travelling waves with dashed curves (red dashed online). We observe that the 
presence of the Hilbert transform increases the value of z/ for which instability arises [70], 
but it does not change the shape of the bifurcation diagram. This is because the Hilbert 
transform term acts as a negative diffusion, see Equation (1.6), and therefore its presence 
acts to shift the bifurcation diagram to higher u, i.e. lower a — 4/z/ as seen in the figure. 
We emphasise that the bifurcation diagrams in Figure 1 are not complete and we expect 
additional unstable branches in analogy with known results for the KS equation [39]. 
This is not a restriction here, since we are interested in demonstrating the stabilisation 
of unstable steady or travelling wave solutions, rather than the stabilisation of all such 
branches. For the branches computed here, we analysed their stability numerically by 
adding a small perturbation to the initial condition (about 10% or smaller of the am¬ 
plitude of the steady state solution) and studied the time evolution to ensure that we 
identified unstable steady solutions to be stabilised using linear feedback controls. 

3.2 Time dependent simulations and feedback control 

We used a Galerkin truncation [67] for the spatial discretisation of the PDE, with the 
number of modes varying between 32, 64 and 128 depending on the number of unstable 
modes. Time integration is carried out using second order implicit-explicit backward 
differentiation formulae (BDF) schemes [1, 2], 

To construct the matrix K necessary for the stabilisation of the steady states, we 
used Matlab 7 s command place. Given the matrices A and F>, we sought a matrix K such 
that the eigenvalues of the matrix A + BK were: 

• —1 if it is the eigenvalue corresponding to the constant eigenfunction -A=. 

\ ZlX 

• ±A if A is an eigenvalue of A with negative/positive real part. 

• — lOhA instead of —A if 6 > 0. We do this because the amplitude of the solutions 
grows with S - [38], so we need to account for this when building the controls. 

We begin by presenting numerical results in the absence of electric fields and dispersion 
(/x = 0, 5 — 0) and for two values of v — 0.2 and v — 0.4 (note that the number of 
unstable eigenvalues is 21 + 1 where l — [z/ -1 / 2 ], where [•] denotes the integer part. The 
number of controls used is 5 and 3, respectively, i.e. equal to 21 + 1; these are placed 
equidistantly and the initial condition is 

1 1 5 

uq(x) — _ H—— V (sin (nx) + cos (nx)). 

V2n V 77 i 
n— 1 

The results are presented in Figure 5 and clearly show that the system is controlled to 
the zero solution long before the final computed time of t = 5; our results are also in 
good agreement with those in [12]. 
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(a) v — 0.2 (b) v — 0.4 


Figure 5: Spatiotemporal evolution showing stabilisation to the zero solution of the KS 
equation for (a) v = 0.2 (a = 20), and (b) v — 0.4 (a = 10). 


Results analogous to those presented in Figure 5 were found regarding the stabilisation 
of the zero solution to the KS equation in the presence of an electric held. In what follows 
we use the following initial condition unless stated otherwise: 

uo(x) = ~^= (sin(x) + cos(x)). (3.7) 


Note that the number of unstable modes is 2/ + 1 where l 


47 

27 


: see Proposition 


1- and this is the number of controls used in the numerical experiments. The numerical 
results for v = 0.2 and fi = 0.5 with 5 equidistant controls are shown in Figure 6, where 
we again clearly observe stabilisation to the zero solution. 



Figure 6: Spatiotemporal evolution showing stabilisation to the zero solution of the KS 
equation in the presence of an electric held with p, = 0.5 and v = 0.2 (a = 20). 

Having shown the stabilisation of zero states for relatively small values of v. we turn 
next to the stabilisation of nontrivial steady states of the generalised KS equation (1.4), 
in the absence of dispersion. We illustrate the feasibility of our control methodology 
for two typical cases that yield unstable steady states as computed in the bifurcation 
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diagram of Figure 1. In the first case we use v — 0.1115, // = 0, and in the second 
v — 0.35, n = 0.3. In both cases we used 2/ +1 equidistant controls, i.e., the same as the 
number of unstable eigenvalues of the system. The results of our numerical experiments 
are presented in Figures 8 and 9, respectively. When v — 0.1115, [i — 0, i.e. a ~ 35.87, 



Figure 7: Zoom in of Panel 1(a) with v £ [0.1,1]. Branches are labelled as used in 
Tables 4-6: Branch 1 unimodal steady states; branch 2 - bimodal steady states; branch 3 
- trimodal steady states; branch 4 - tetramodal steady states. The cross and open circle 
symbols indicate the steady states (stable and unstable, respectively) that are shown in 
the Figure 8. 

both stable and unstable steady states coexist and the solution of the PDE with a given 
initial condition, e.g. (3.7), evolves to the most attracting stable state. This is shown 
in figure 8(a) where it is seen that the solution evolves to a stable bimodal steady state, 
marked with a circle in Figure 7. We are interested in using feedback control to stabilise 
one of the coexisting unstable steady states, and the results of achieving this are presented 
in Figures 8(b)-(c); panel 8(b) shows the evolution of the initial condition (3.7) using 
2/ + 1 = 5 equidistant controls and stabilisation of the steady state marked with a + in 
Figure 7 is achieved relatively quickly after approximately 2 time units. The evolution 
of the amplitudes of the 5 applied controls is shown in Figure 8(c), and we see that the 
required energy tends to values very close to zero as time evolves. Note that the control 
amplitudes remain small and close to zero once the unstable controlled state is reached, 
but they cannot be identically zero due to the unstable nature of the controlled solution. 
Figure 9 shows the results for v — 0.35, fi — 0.3. The solution we choose to stabilise at 
these values is an unstable bimodal steady state and Figure 9 shows how it is stabilised 
using 21 + 1 = 5 controls. 

Our last task is to stabilise travelling wave solutions of the equation with and without 
dispersion and electric held. Figure 10 illustrates the stabilisation of three different 
travelling wave solutions to the KS equation (1.4) with no dispersion or electric held 
(5 — fi = 0) and for a small value of v {y = 0.01) which corresponds to a very large 
domain (L — 20n « 62) that enables the existence of single pulse travelling waves as well 
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(a) Uncontrolled Solution (b) Controlled Solution 



Figure 8: Control of non-uniform solutions of the KS equation for v — 0.1115; panel (a) 
spatiotemporal evolution without controls (solution belongs to branch 1 of the bifurca¬ 
tion diagram in Figure la); panel (b) controlled to the steady state in branch 4 of the 
bifurcation diagram in Figure la; panel (c) evolution of the amplitude of the 5 applied 
controls. 



Figure 9: Spatiotemporal evolution of the stabilised steady state of the Kuramoto- 
Sivashinsky equation for v = 0.35 (a ~ 11.43), fi — 0.3 . 


as two- or three-pulse bound states. However, due to the small value of v. when solving 
the PDE the initial condition evolves to a solution that exhibits the spatiotemporal 
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chaotic behaviour that is characteristic of this equation. Panel 10(a) shows this chaotic 
behaviour while Panels (b)-(d) show the evolution of the controlled solution to 1, 2 and 
3 pulses, respectively. We used m — 21 equidistant controls in each case. 




Figure 10: Solution to the KS equation for v — 0.01 (10a) with no controls, and controlled 
to (10b) one solitary pulse, (10c) two solitary pulses and (10d) three solitary pulses. 

The control algorithms presented and analysed here also apply to the gKS equation 
with 5 > 0. The findings are similar and for brevity we do not present them here. 
Detailed results and animations can be found in [25] and its supplemental material. 

4 Optimal control for the generalised Kuramoto-Sivashinsky equation 

In practical applications it is important to apply controls that minimise the cost asso¬ 
ciated with their use, and this leads to considerations of an optimal control problem 
based on some measure of the energy cost of the controls. In what follows we consider 
the energy of the controls given by their T 2 -norm. Since the controls decay to almost 
zero relatively fast in time (see Figure 9, for instance), we expect that minimising their 
L 2 -norm should decrease their amplitude. 

The objective, then, is to achieve control of the zero solution or other unstable steady- 
state or travelling wave solutions of the gKS, and to do this while spending the least 
energy possible. To that end, we consider a cost functional that includes the distance 
between the solution and the desired state, as well as the L 2 -norm of the controls used. 
Different distance norms ||it — u|| can be used and in our computations we employ the 
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L 2 , H 1 or H 2 norms. The reason we consider different norms is that the solutions are 
expected to belong to H 2 ( 0, 2n), and we wish to analyse the effects of the regularity and 
the oscillations of the solutions on the cost functional. 

In flow control it is possible to use point actuator functions [4, 11, 12, 21], implying 
that we can take bi{x ) = 5(x — Xi). Assuming that the cost of placing a control at 
Xi is the same for all actuator positions x G (0, 27 t), it makes sense to seek a solution 
that minimises the norms of the control functions fi(t). Since the delta functions in the 
feedback controls are not L 2 functions, the standard results of constrained optimisation 
for PDEs [44, 68 ] do not apply. Because of this hurdle we will first prove existence of 
optimal controls for the case of general controls, f(x,t ) G T 2 (0, T; T 2 (0, 27 t)), i.e. mean 
zero spatially periodic controls in T 2 (0, 27t) that are also L 2 functions of time, and focus 
on the case of feedback controls where we can apply standard optimisation techniques. 

We consider cost functionals of the form 

rj~\ rj~i JYl 

c ( u , f ) = 7)[ IN-,*) - u\\ 2 dt+ i||u(-,T)- m|| 2 +! f ^fi{t) 2 dt (4.1) 
do do i=1 

where u is the desired steady state, F = [f\ • • • f m \, and the norm (e.g. L 2 , H 1 or H 2 ) 
is left unspecified. The choice of the parameter 7 depends on how large we are willing 
to allow the norm of the controls to become: if we need to maintain a small norm of the 
controls while allowing the solution to be considerably different from the steady state, 
then we use 7 > 1. If, on the other hand, we have a very large amount of energy to spend 
on the controls and want the solution to be as close as possible to the desired steady state, 
then we choose 7 << 1 so that the weight of the controls does not influence significantly 
the value of the cost functional. The terminal time term ^||u(-, T ) — u || 2 is introduced to 
provide us with a condition for the final-value problem obtained when solving the adjoint 
equation for the optimisation problem. Before proceeding to the optimisation problem 
it is useful to introduce the following definitions. 

Definition 1. The space of admissible contols, F ad , is a bounded convex subset of 
L 2 (0,T;L 2 (0,27t)). 

Definition 2. A control f* G F ad is said to be optimal, and u* — u(f*) is the associated 
optimal state, if C(u(f*),u, /*) <C(u(f),u,f ) V/ G F ad . 

Our numerical experiments presented in Section 4.2 suggest that, given an initial 
condition and a desired steady state, there exists at least one optimal placement of the 
control actuators for every value of v and /x. However, here we prove existence of an 
optimal control in the case of an open-loop control using the quadratic cost functional 

1 1 '■y f 27T 

C ( U J) = 2 J o ll u (’’ t ) “ ^\\h dt + 2 lK'’ T ) + l J o J o f(x,t) 2 dxdt. (4.2) 

The point actuated controls in the form of delta functions are not in L 2 and hence an 
analogous proof in this case requires distribution theory which is beyond the scope of the 


22 



present study. The optimisation problem is: 


minimise C (it, /) 

(4.3a) 

subject to u t + vu xxxx + u xx + uu x = f(x, t), 

(4.3b) 

u(x, 0) = ito(x) € Hp( 0, 27t), 

(4.3c) 

(Fu d^u. N . ^ ^ „ _ 

— j (x + 2*) = ^(x), j = 0,1, 2,3, 

(4.3d) 

/ e F ad . 

(4.3e) 


The main result of this section is the following. 

Theorem 1. Assume that F a d c T 2 (0, T; T 2 (0, 2n)). Then (4.3) has at least one optimal 
control /* with associated optimal state u*. 

Remark 6 . Since the Hilbert transform and third derivative terms are linear functionals 
of u, Theorem 1 can be easily generalised to the case when //, <5 > 0. 

Remark 7. The presence of the Burgers nonlinearity in the PDE makes the optimisation 
problem no longer convex. Consequently, we do not expect the solution of the optimal 
control problem to be unique. 

The nonlinearity in our problem, defined by A/(it) — uu x , is twice Frechet differen¬ 
tiable with respect to u but is neither an increasing functional of u, nor is it globally 
Lipschitz continuous. Furthermore, it depends explicitly on the derivative u x . Conse¬ 
quently, the well developed theory of optimal control for systems of reaction-diffusion 
equations [44, 68 ] does not apply to our problem; see Equation (4.3). 

Proof of Theorem 1. Let X — iL 1 (0, T; H 2 ( 0, 2n)) x F a d and e(-, •) be a functional defined 
in X by 

Ut + nuxxxx T u X x T uu x / 
it(-, 0) — uq{x) 

Our optimisation problem is that of minimising the cost functional C subject to e(tt, f ) = 

0, periodic boundary conditions and (it,/) £ X;; see Equation (4.3). 

Let (it, /) £ X satisfy e(it, /) = 0. Since C is a function of the sum of the norms of u 
and /, it is clear that C is nonnegative and 

C{u, /)->■ oo for ||(it, f)\\x oo. (4.5) 

Therefore, there exists a constant c > 0 such that c — inf e ( M j ) =0 C(it, /) = lim^oo C(u n , f n ), 
where (it n ,/ n ) is a minimising sequence in X , which exists due to the reflexivity of L 2 . 
From Equation (4.5) we can conclude that {(it n , / n )} neN is bounded, and therefore there 
exists (it*,/*) € X such that (it n ,/ n ) —^ (it*,/*) for n —> oo. This means that all 
the linear functionals of it n and f n , and in particular their derivatives, also converge 
weakly to the same functionals of it* and /* in the appropriate space. Hence, we only 
need to prove the convergence of the nonlinearity. 


e(u,f) 


(4.4) 
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Following an argument similar to that in [73] for the Burgers equation, we notice 
that since for every t G [0 ,T] we have u*(-,t) G H 2 ( 0,27t), then € C'([0, 27 t]) and 

therefore if ip G X, (u*(p)(-,t) G T 2 ([0,27 t]). Hence, w*<p G i7 1 (0, T; L 2 (0, 27t)) and 

pT p2tt 

/ / (u n -u*) x u* V dxdt— w0, (4.6) 

do do 

Finally, from the estimates (2.1) we know that ||u™|| is bounded, and since H 2 (Q) is 
compactly embedded in T 2 (fl), we deduce that 

pT p2tt 

/ / (w n - cte dt < \\u n - U*|| IKII ||<^|U°o( 0 , 27 r) - 0 , 

do do 

e H?,(0,2tt). (4.7) 

Hence, by adding and subtracting appropriate terms we have 

pT p2it p2tt 

/ / (u n u™ — u*u*)(p dx dt = / {u n — u*)u™(p + (u n — u*) x u*(p dx —0, 

do do do 

v*7 e ffp(0, 2ff), (4.8) 

and therefore the nonlinearity u n u™ is weakly convergent to u*u* in X. Now, noticing 
that u* and u* are continuous in [0, 27 t] x [0, T], we observe that u* satishes the periodic 
boundary conditions and the initial condition. If we now consider ip G X satisfying 
ip(x,T) = 0, and use the weak convergence of the derivatives of u and equation (4.8), 
we conclude that (u*, /*) is a weak solution of the state equation. The optimality of the 
pair (u*, /*) follows from the weak lower semi-continuity of C (cf. proof of Theorem 4.15 
in [ 68 ]). □ 

4.1 Algorithm and Numerical Experiments 

We note that the dependence of the cost functional on the positions X{, i = 1,..., m is 
in the controls, since the matrix K necessary to define them depends on the positions 
chosen. However, when defining the Lagrangian we will assume that only the functions 
bi(x) depend on X{, and treat the controls /)(t) as if they were independent of the positions 
Xi. Under this assumption we are able to obtain very satisfactory results, as evidenced 
by the results presented in the tables below. We begin by introducing the Lagrangian 


C (u,p, [xi,x 2 , ■ ■ .,Xm] T ) = 


| u(-,t) — u\\ 2 dt + -||u(-,T) — u\ 


dt 


(u t + VU 2 


Uxx + uu x ) p{x , t) dx dt 


'o Jo 

pT p2tt m 


5(x — Xi)ifi(t)p(x, t ) dx dt. 
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Integrating by parts in space and time and computing the Frechet derivative with respect 
to u (with test functions h(x,t ) satisfying h(x, 0) = 0), we obtain the adjoint equation 


Pt “1“ VPxxxX “1“ Pxx UPx 
p(x,T) 

0( x + 27r ) 


YjT= 1 K x ~ Xi)I<i.Z p + U-U, 
u(x. T ), 

Sw. 


(4.10) 


where x E [0, 2n] and t E [0,T]. This PDE is backwards in time but is well-posed 
since it is a final value problem. To solve it, we obtain the discretised ODE system 
—z p — Az p + G ad j(z p , z u ) + z u — z u , where the elements of G ad ^(z p , z u ) are given by 


1n,adj 


9n,adj 


2 ^ j+k=n 


Y H u jPk - u jPk ) 


2 ^ . 


Y ( k ( u jPk + u ]pd - j( u kPj + u kP c j )), 


j—k=n 


Y k ( 1J pPk + u jPk) + 2^/n ( k ( u jPk - u jPk) + j(u C kPj - u kPj)) 


j+k=n 


j—k=n 


and we have used the Fourier series representation p(x,t) — -^ 7 = + Pn(t) S 11 ^k x ' > + 

F^°° cos(nx) 

2 ^n=lPn\ L ) ■ 

Differentiating with respect to the positions of the control actuators, we also obtain 
a descent direction using the variational inequality, or first variation 

rT 

/ [fl(t)Px(xi,t) ■ ■ ■ fm(t)Px(x m ,t)} T ’ (x ~ x) dt > 0, Vx = [x 1 • • • X m ] T , (4.11) 
Jo 

where x — [aq,... , x m ] are the optimal positions. To proceed with the optimisation, we 
will use a gradient descent method, see [68, Sec. 5.9], and consider F ar ] = (0,27r) m . The 
algorithm is as follows. 

Note that as mentioned earlier we consider the following three different cost function¬ 
als: 


Ci (u,u,f) = - 


|w(-, t) — u || 2 dt 


-||w(-,T) -u 


r: 11 2 


+ li Y , ii/^(^)Hl 2 (o,t) ( 4 - 12 ) 


i =1 


C 2 (u, u, f) = i/ 0 T (||u(-,t) - u || 2 + ||w x (-,t) - w x || 2 ) dt 

m 

+ 2 ll^( i )ll^ 2 (0,T) 

i= 1 

C 3 (w, W, y) — 2 Jo (ll^(’> ^11 ~F ll’^'x(’) ^x|| T ||Wxx(‘j ^xx|| ) dt 

+ 2 — ^l| 2 ll' U ic(’;^) — Wx|| 2 + |Kx(-,r) — W^ll 2 ) 

m 

2 = 1 


(4.13) 


(4.14) 
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Algorithm for optimal control for the KS equation 

Given u, /x, 7 , T, uo(x), A, x°, B°, u, compute the matrix K°. 
while F(current iteration) < C (previous iteration) do 

1. Solve the state equation to obtain u k ~ x and compute F(xx fc_1 , xx, F{x k ~ x ))\ 

2 . Solve the adjoint equation to obtain p k ~ 1 ; 

3. Define P, = p k ~ 1 (xp 1 ,t) ■ ■ ■ , 

Pk- 1 = [ K k ~ 1 (z v,k ~ 1 — z™)P x dt and hk = —Pk-i; 

Jo 

4. Find s — min s> o { C(u(x k ~ 1 + shk ), xx, F{x k ~ 1 + shk ))}; 

5. Project x k ~ x + shk into (0,27r) m , obtaining x k ; 

6 . Compute the matrix B k \ 

7. Compute the matrix K k with MatlaF s command place. 

end 

Figure 11: Algorithm for optimal control of the KS equation. 


26 




4.2 Numerical Experiments 

Computations were carried out using the algorithmpresented in Figure ?? for various 
values of v and /x. The number of controls used was equal to the number of unstable 
eigenvalues, and m equidistant points were used as an initial guess for the position of the 
controls. In all the computations the initial condition is uq { x ) — S1 ^' > + cc ^£y , and the 
final time is T — 10. Numerical results are presented in Tables 1-3 for the stabilisation 
of the zero solution of the KS equation, and in Tables 4-6 for the stabilisation on non¬ 
trivial unstable steady states as computed in the bifurcation diagram of Figure la. Each 
entry in Tables 1-3 contains the value of v. the value of the cost functional (C\. C 2 and 
C 3 for tables 1, 2 and 3, respectively), the cost of the controls X^=i ll/*(^)H l 2 (o,t)> ^he 
number of iterations required to obtain an optimal state, and in the last column the 
spatial distribution of the controls over the domain [ 0 , 2n] - a heavy dot is placed where 
a control acts. Tables 4-6 are presented in an analogous manner, with the difference that 
the first column provides information on the unstable solution that is being controlled, 
and in particular the branch on Figure la where the solution was taken from is stated 
along with the value of v. As the results indicate, several distinct unstable solutions at 
a given value of v are controlled (e.g. for z/ = 0 . 1 , three solutions are stabilised coming 
from branches 1, 3 and 4 respectively). 


V 

Cost Ci 

Cost of Controls 

Iterations Optimal Positions 

0.9 

8.9647 

0.5592 

q L.L -x-x-x-x_!-•-x- 

O 0 1 2 3 4 5 6 

0.8 

6.5012 

1.1274 

<3 E-•-*-*-*-i-•-*-x— 

U 0 1 2 3 4 5 6 

0.7 

5.5760 

1.7204 

A E - • -X-X-X-±1- • -X-X- 

L ± 0 1 2 3 4 5 6 

0.6 

5.2803 

2.3018 

3 0 12 3 4 5 6 

0.5 

5.2230 

2.8204 

r E-•-x-x-x-• •-x- 

O 0 1 2 3 4 5 6 

0.4 

5.9339 

3.8404 

3 0 1 2 3 4 5 6 

0.3 

6.2152 

3.9813 

O E-•- X-X-X-• -• -x_ 

Z 0 1 2 3 4 5 6 

0.2 

6.3127 

4.5261 

9 E-5_i- 1 - 1 -*-1-x_?-*-'-x— 

Z 0 1 2 3 4 5 6 

0.1 

7.1652 

5.5759 

Z 0 1 2 3 4 5 6 


Table 1: Optimal positions and value of the cost functional considered in the L 2 -norm 
for different values of v when stabilising the zero solution to the KS equation. 

As expected we observe that the value of the cost functionals C 1 , C 2 , C 3 given by 
(4.12)-(4.14), increases as v decreases. Furthermore, the value of the cost functional also 
increases as we increase the desired regularity of the solution from L 2 to H 1 to H 2 — norms. 
Comparing the results and in particular the positions of the optimal controls for the three 
different cost functionals in Tables 1-3, we can conclude that for the stabilisation of the 
zero steady states, the optimal control problem is more robust (in the sense that the 
optimal positions of the controls do not change much as v is reduced) when the L 2 cost 
functional C\ is used. 

Turning now to the results of Tables 4-6 that deal with the stabilisation of unstable 
nonuniform steady states, we observe once again that there is an increase in the cost 
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V 

Cost C 2 

Cost of Controls 

Iterations 

Optimal Positions 

0.9 

17.6354 

0.9698 

3 

E_1_1_L?_1_ * _1_1_• 

0 1 2 3 4 5 6 

0.8 

12.2553 

1.5486 

6 

0 1 2 3 4 5 6 

0.7 

10.8270 

2.0016 

2 

0 1 2 3 4 5 6 

0.6 

10.2340 

4.0181 

6 

0 1 2 3 4 5 6 

0.5 

10.1517 

4.5407 

4 

0 1 2 3 4 5 6 

0.4 

9.3345 

4.2298 

2 

0 1 2 3 4 5 6 

0.3 

10.9671 

4.8110 

2 

0 1 2 3 4 5 6 

0.2 

10.7154 

5.5308 

3 

0 1 2 3 4 5 6 

0.1 

9.7088 

5.6463 

1 

0 1 2 3 4 5 6 


Table 2: Optimal positions and value of the cost functional considered in the iiO-norm 
for different values of v when stabilising the zero solution to the KS equation. 


V 

Cost C 3 

Cost of Controls 

Iterations 

Optimal Positions 

0.9 

27.2098 

1.1313 

3 

0 1 2 3 4 5 6 

0.8 

19.3431 

2.3490 

5 

E_1_e_1_1_ • _1_1_• 

0 1 2 3 4 5 6 

0.7 

15.8522 

2.5815 

4 

0 1 2 3 4 5 6 

0.6 

14.0865 

3.3384 

5 

0 1 2 3 4 5 6 

0.5 

17.0462 

6.4166 

1 

0 1 2 3 4 5 6 

0.4 

20.7720 

8.3217 

1 

0 1 2 3 4 5 6 

0.3 

14.4393 

5.3865 

3 

0 1 2 3 4 5 6 

0.2 

21.5856 

6.1456 

1 

0 1 2 3 4 5 6 

0.1 

21.1636 

5.6463 

1 

0 1 2 3 4 5 6 


Table 3: Optimal positions and value of the cost functional considered in the H 2 -norm 
for different values of v when stabilising the zero solution to the KS equation. 


V 

Cost C\ 

Cost of Controls 

Iterations 

Optimal Positions 

0.3 

14.3164 

9.7419 

5 

0 1 2 3 4 5 6 

0 . 2 ,br.l 

30.0588 

21.8691 

1 

0 1 2 3 4 5 6 

0.2,br.3 

24.0520 

16.2832 

2 

0 1 2 3 4 5 6 

0 . 1 ,br.l 

28.5591 

32.9859 

2 

0 1 2 3 4 5 6 

0.1,br.3 

37.8902 

32.4264 

1 

0 1 2 3 4 5 6 

0.1,br.4 

62.3916 

51.6820 

4 

0 1 2 3 4 5 6 


Table 4: Optimal positions and value of the cost functional considered in the L 2 -norm 
for different values of v when stabilising some of the nontrivial steady states from the 
bifurcation diagram la. 


functionals as v decreases. We also observe that in this case (and in contrast to the 
stabilisation of the zero solution) the higher order norms give optimal controls that are 
more robust, with respect to changing v. in comparison to utilising the L 2 cost functional. 
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V 

Cost C 2 

Cost of Controls 

Iterations 

Optimal Positions 

0.3 

24.6922 

9.6582 

3 

E_•_._._•_._ 

0 1 2 3 4 5 6 

0 . 2 ,br.l 

53.7672 

17.9929 

2 

0 1 2 3 4 5 6 

0.2,br.3 

52.2630 

15.4089 

3 

0 1 2 3 4 5 6 

0 . 1 ,brl 

87.1636 

35.2169 

1 

0 1 2 3 4 5 6 

0.1,br.3 

83.9787 

33.7581 

2 

0 1 2 3 4 5 6 

0.1,br.4 

171.6040 

62.1381 

2 

0 1 2 3 4 5 6 


Table 5: Optimal positions and value of the cost functional considered in the iiO-norm 
for different values of v when stabilising some of the nontrivial steady states from the 
bifurcation diagram la. 


V 

Cost C 3 

Cost of Controls 

Iterations 

Optimal Positions 

0.3 

54.5441 

13.8436 

2 

0 1 2 3 4 5 6 

0 . 2 ,br.l 

262.7363 

21.1679 

3 

0 1 2 3 4 5 6 

0.2,br.3 

266.0515 

32.4603 

3 

0 1 2 3 4 5 6 

0 . 1 ,br.l 

702.4697 

35.2169 

1 

0 1 2 3 4 5 6 

0.1,br.3 

745.6007 

32.4264 

1 

0 1 2 3 4 5 6 

0.1,br.4 

1384.6689 

63.1272 

2 

0 1 2 3 4 5 6 


Table 6 : Optimal positions and value of the cost functional considered inthe H 2 -norm 
for different values of v when stabilising some of the nontrivial steady states from the 
bifurcation diagram la. 


Similar numerical experiments were performed for the optimal control problem for 
the KS equation in the presence of an electric field, p > 0. The results are quite similar 
to the ones already presented in this section and we omit providing additional graphs 
and tables. However, we summarise the conclusions drawn from the nonzero electric 
held numerical experiments as follows: with few exceptions, an increase in the intensity 
of the electric held parameter p increases the cost of the controls. In addition, it is 
found that the optimal controls for stabilising zero steady states are more robust, with 
respect to changes in p, when using the L 2 cost functional. Similarly, when stabilising 
nontrivial steady states, more robust optimal positions for the controls arise when the 
H 1 and H 2 cost functionals are used. Both of these hndings are analogous to those for 
the non-electrihed control problem. 

A more detailed comparison of the energy required to control different solutions using 
equidistant actuators or optimally computed positions as described above, is provided in 
Figures 12 and 13. Figure 12 shows the stabilisation of a nonuniform steady state for 
the KS equation, v = 0.3 and p — 0, while Figure 13 shows analogous results but for 
the electrified problem with parameters v — 0.5 and p = 0.4 (in both cases dispersion 
is absent, <5 = 0). Panel (a) shows the spatiotemporal evolution to the desired state 
in the presence of controls, while panels (b) and (c) depict the evolution of the control 
amplitudes (there are three controls in each case) for equidistant or optimally positioned 
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(a) Controlled solution 




Figure 12: Controlled steady state of the KS equation for v = 0.3 (12a) and controls 
applied: (12b) equidistant and (12c) optimal. 


actuators, respectively. The results show that the amplitudes of optimally placed controls 
decay to zero faster than those of the equidistantly placed ones. 

4.2.1 Optimal control of travelling waves 

We also performed similar numerical experiments to find the optimal position of the 
control actuators when stabilising travelling waves. We found that in most cases, we 
cannot do better than equidistant controls. 

We believe that this is due to the following reasons. First, the length of the domain 
needed for the existence of (unstable) travelling waves is long, and therefore the number 
of unstable modes (and hence of the number of controls) is large (e.g., in the example 
of Fig. 10 we are using m = 21 controls); thus, shifting the position of the controls in a 
relatively large domain should not have a big effect on their amplitudes. Second, solitary 
pulses on long domains necessarily have large flat regions which are susceptible to linear 
instabilities leading to the nonlinear wavy perturbations seen in panels (b) and (c) in 
Figure 2. lit is interesting to note that there are 10 wavy structures corresponding to 
the number of linearly unstable modes; thus, we expect optimality when the controls are 
approximately equally spaced thus guaranteeing one control under each wavy structure. 
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(a) Controlled solution 




Figure 13: Controlled steady state of the KS equation for v — 0.5, fi — 0.4 (12a) and 
controls applied: (13b) equidistant and (13c) optimal. 


Shifting the controls can introduce instability and nonlinear growth to a different state. 

5 Feedback and optimal control for coupled Kuramoto-Sivashinsky equa¬ 
tions 

In Sections 2-4 we studied analytically and computationally the feedback control and 
optimal control problems for the generalised Kuramoto-Sivashinsky equation. In appli¬ 
cations, systems of equations emerge with two or more nonlinear coupled PDEs, and 
this section is concerned with the control of such systems. As an example we refer to 
the system of two coupled KS equations that arises in the weakly nonlinear asymptotic 
analysis of a three layer flow of immiscible viscous fluids stratified in a channel and driven 
by gravity and/or a stream wise pressure gradient - see [49]. The fully coupled system is 
a challenging PDE problem and questions such as existence and uniqueness of solutions, 
steady states and bifurcation theory are still poorly understood. Such problems are cur¬ 
rently under investigations and our findings will be reported elsewhere. In this Section we 
consider the problem of feedback and of optimal control for a system of KS equations that 
are coupled only through the second derivatives. Such a coupling is special but can arise 
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iii the application of three-layer flow. More generally, the nonlinearities are also coupled 
and in fact the nonlinear flux functions can generically have real or complex eigenvalues 
implying hyperbolic elliptic transitions thus complicating the analysis significantly; see 
[48] for a detailed study of such effects. 

In what follows we consider the following coupled system of Kuramoto-Sivashinsky 
equations 

\ ^l,t — VU1,XXXX ^1,XX W Ul fX <^1^2, XX /r -i \ 

\ ^2,t — VU2,xxxx ^2,xx ^2^2,x ^2^1,xx- 

We consider the equations in the interval (0, 27 t) with periodic boundary conditions and 
initial conditions m(%, 0) = mo(%) and U2(x,0) — U 2 q{x) and wio, W 20 £ H 2 ( 0, 27t). 

We can prove, using the background flow method [17, 46, 71] that the solutions to 
the system are bounded: 

Proposition 5. Assume that mo, U 20 E 77^(0, 2-7t). Then there exists a constant C — 
C(u,a 1 , 02 ) such that 

||wi||l 2 + ||u2||l 2 — C . (5-2) 

Similar bounds can be obtained for the 77 1 - and H 2 -norm, and we can therefore conclude 
that the solutions to system (5.1) are in L°°(0, 2-7t). 

The proof is not included here due to space constraints and will appear elsewhere. A 
similar result can be proved for the case when the coupling comes only through the fourth 
order derivatives (i.e. there is a non-diagonal negative definite fourth order viscosity 
matrix), and we believe that similar results can also be proved for the fully coupled 
system with non-diagonal second as well fourth order viscosity matrices. 


5.1 Feedback control for the coupled KS equations 


Since Equations (5.1) are coupled linearly, analogous results to the ones presented earlier 
for the single KS are obtained. First, we can prove that it is possible to stabilise any 
steady state solution (either the zero solution or any nontrivial steady state) for this 
system. We proceed in the same way as for the single KS equation and write the con¬ 
trolled system (for completeness we have also included coupling through the fourth order 
derivatives) 


r — 

1 w 2 ,t — 

Defining 

U (x, t ) = 


VUi^xxxx ^l,xx UiUi,x ®lU2,xx + $(x Xj 1 )fji (^) ’ 

~UU 2 ,xxxx ^2 ,xx 'U j 2^ j 2,x C^2^1,xx T — i d(x Xj, 2 ) f j 2 (^) • 


1 

/ ^ s 

jT 

s 

1_ 

00 

_ v 

’ uf n (t) ' 

1 

s 

_1 

n =1 

. u s 2n (t) _ 


00 

sin (nx) + 

n =0 


n (t) 
U°2 n (t) 


cos (nx), 


(5.3) 


(5.4) 


and taking the inner product with the functions -^==, sll ^T r ' ) and cos (fh x l yields the fol¬ 
lowing infinite system of ODEs 


= (- un4 + n2 ) U in + a i n2uS jn + 9in + Y,Ti=l b iinfji( t ) n=l,..., 00 , 


= (-^ 4 + n 2 ) u c in + a in 2 u ( } n + g c n + TX=i h lJn 0) n = 0,...,o o, 


(5.5) 
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where i, j = 1,2 , i ^ j, and the functions b and g are dehned as in the single KS case. 
We truncate the system at N modes and define 


£ 


u 


— [u\ o uh u\i ■ ■ ■ U S 1N u\ N U20 U21 U21 ■ ■ ■ U2N u\n] T , 
G — [0 g\ 1 Qii • • • gf N g\ N 0 g^i g%i • • • 9 2 n 9 2 n] T > 


F = [fll(t) fu(t) ••• flm(t) f2l(t) f22(t) ••• f 2 m(t)] T . 


Next we write 


A = 


A() A\ 

A 0 ’ 


B = 


B 1 

B 2 


where Aq — diag(0, — v + 1, — v + 1, • • • , — vrA + n 2 , — un 4 + n 2 , • • •), Ai 
oti, -Pi + Oii, ■ ■ ■ , - Pin 4 + OLin 2 , - Pin 4 + a*n 2 , • • •) and 


diag(0, -pi + 


B x = 


b 

b 

b 


c 

1*0 

h c 

°2i0 

• • b c 

u rrii0 

s 

1*1 

b s 

0 2,1 

■ ■ b s 

y mil 

c 

1*1 

b c 

°2j 1 

■ ■ b c 

^TOjl 


for i — 1,2. Hence the inhnite system of ODEs can be written as 

z u = Az u + G + BF. (5.6) 

We can prove a result similar to Proposition 1. 

be a nontrivial (unstable) steady state solution of (5.1), 

and let l — l\ + l 2 be the number of unstable eigenvalues of the linearised system, i.e. 
l\ < 1 + W * lQ;2 < (^ _|_ if 2 anc [ i 2 < 1 VAPAPl <; (/ 2 4 - 1 ) 2 . if m = 2(1 + 1 ) and there exists 
a matrix K such that all of the eigenvalues of the matrix A + BK have negative real part, 
then the state feedback controls 

[flip) fl2(t) ••• flm(t) f2l(t) f 2 2(t) ••• f2m(t)] T = F = K(z U ~ Z°), (5.7) 

stabilise this nontrivial steady state solution of system (5.1). 

The proof of this result follows the same argument as the proof of Proposition 1. 

We present in Figures 14 and 15 the numerical results of the stabilisation of the 
zero solution and a steady state solution, respectively, for system (5.1) with v — 0.5, 
a\ — 0.8 and a 2 — 0.5. We used m — 4 equidistant controls to control each solution, 
corresponding physically to applying 4 controls in each wall. Upper panels correspond 
to the uncontrolled solution, and lower panels correspond to the stabilised solution. We 
clearly observe in both figures the stabilisation of the desired steady state. 


Proposition 6. Let U — 


U\ 

U 2 
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Figure 14: Uncontrolled solution (v{, i — 1,2) and controlled zero solution (tq, i = 1,2) 
of the system of coupled KS equations for v — 0.5, aq = 0.8 and a 2 = 0.5. 



Figure 15: Uncontrolled solution (vi, i — 1,2) and stabilised steady state solution (iq, 
« = 1,2) of the system of coupled KS equations for v = 0.5, a\ = 0.8 and 0:2 = 0.5. 
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5.2 Optimal control of the system of coupled KS equations 


We consider next the problem of controlling an arbitrary steady state U = 
optimal way. We introduce the cost functional 

"T 


1 


C{U 7 F)= -J (IM-,*) - ui \\ 2 L 2 + \\u 2 (-,t) - U 2 WI 2 ) dt 
+ i(||i 11 (,T)-S I ||l 2 + |M,T)-S 2 ||l 2 ) 

+ 'i Jo dt - 


U\ 

U 2 


m an 


(5.8) 


The optimisation problem that we have to solve takes the form 


minimise C (U, F) (5.9a) 

subject to ui t t T T u\, X x T uiui,x T ey.\U 2 , X x — /1 ((£?£)? (5.9b) 

^2 ,t T L / U2,xxxx T U2,xx T ^2^2,x T 0 ) 2^1 ,xx = f2{x,t), (5.9c) 

Ui(x, 0) = uo t i(x), i — 1,2, (5.9d) 

F^-(x + 2it) = F^-(x), j = 0,1,2,3, * = 1,2, (5.9e) 

fieF ad , * = 1,2. (5.9f) 


Here, uq ^ G HJ(0,2tt) and F ac [ is a bounded, closed and convex subset of T 2 ((0, 2n) x 
(0 ,T)). ’ 

We can prove the following theorem. 


Theorem 2. If F cu j c T 2 ((0, T); L 2 (0, 27t)) ; the optimal control problem (5.9a)-(5.9f) has 


at least one optimal control F* — 


ft 

ft 


with associated optimal state U* — 


Sketch of the proof. The proof follows the same steps as that of Theorem 1 for 
the single Kuramoto-Sivashinsky equation. For the coupled system, we need to consider 

a different state space X — (h 1 { 0, T; i? 2 (0, 2^))) x (F a d) 2 , and redefine e(-, •; ■, •): 


e(ui,u 2 ;/i,/ 2 ) 


U\ T T Ul,xx T U\U\^x T eX\U2^xx fl(x,t) 

U 2 ,t T L /r a2,xxxx T U2, xx T U 2 'U 2 j 2 ; T OL2^1, xx f 2 (x,t) 

ui(-,0) - u 0 ,i(x) 
u 2 (-,0) - U 0:2 (x) 


(5.10) 


The rest of the proof follows Theorem 1, but accounting for the fact that for every 
t G [0, T} we have [/*(•, t) G (^Hp{ 0, 27t)^ , and then U*(-, t ) G (C([0, 27t])) 2 and therefore 
if ipi G X, (u*(pi)(- : t) G T 2 ([0, 2ir]) for i = 1,2. 
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Filially, we also need the estimates in Proposition 5 to establish that ||n ^ x ||^2 is 
bounded, and since H 2 is compactly embedded in L 2 , we deduce that 



dx dt < IK 


i,x l 


0, \/tpi G tf 2 (D). 

(5.11) 


6 Conclusions 

In this paper we studied the problem of controlling and stabilising solutions to the gen¬ 
eralised Kuramoto-Sivashinsky equation. We studied both feedback and optimal control 
problems. For the optimal control problem, we proved existence of an optimal control and 
we investigated numerically the problem of optimal actuator placement. By extending 
earlier work by Christohdes et al. [5, 6, 10, 12] we showed rigorously that we can control 
arbitrary nontrivial steady states of the Kuramoto-Sivashinsky equation, including trav¬ 
elling wave solutions, using only a finite number of point actuators. The number of point 
actuators needed is related to the number of unstable modes. We also investigated the 
robustness of the controllers with respect to changing the parameters in the equation. In 
particular, we showed that our proposed control methodology can be used in the presence 
of uncertainty. Our results can be extended to coupled systems of Kuramoto-Sivashinsky 
equations. 

We have not discussed about the practical implementation of the control method¬ 
ologies studied in this paper. For example, we have assumed that complete information 
about the solution of the KS equation is available. This, and other issues related to the 
implementation of the control algorithm are discussed in [66]. 

We considered the case where the entire solution to the KS equation is available to 
us. It is straightforward, however, to apply our results to the case when only a finite 
number of observations is available, using the techniques presented in [5]. For brevity of 
exposition we have not done this for the KS equation. In a forthcoming paper [66] we 
study the feedback control problem when only a finite number of observations is available 
to us for a more complicated PDE, the quasilinear Benney equation arising in falling him 
problems (see [34] and references therein). 

There are several directions in which the results presented in this paper can be ex¬ 
tended. First, the KS equation is a simplified model for thin him hows obtained using 
weakly nonlinear analysis and valid close to criticality [20, 34], We can apply the control 
methodologies studied in this paper to other simplified models that are closer to the 
full 2D Navier-Stokes dynamics, such as the Benney equation and the weighted resid¬ 
ual model. These equations are more complicated since they are quasilinear and they 
can include additional degrees of freedom. It is possible to extend our results to such 
models [66]. Furthermore, the control of unstable travelling waves for the KS equation 
with dispersion can be analysed in detail, leading to an efficient and robust algorithm. 
This problem is studied further in [25]. Finally, our techniques can be extended so that 
they apply to the noisy KS equation. Given that the noise itself can sometimes stabilise 


36 



linearly unstable solutions [54, 55], the interaction between noise and controls can lead 
to very interesting dynamic phenomena. We think that this is a particularly interesting 
direction for further research, since the noisy KS equation is closely related to the noisy 
KPZ equation, which is a universal model for weakly asymmetric processes [28]. 
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